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ABSTRACT 

The possibility is explored that accretion on an intermediate mass black hole contributes to the 
ionisation of the interstellar medium of the Compact Blue Dwarf galaxy MRK 996. Chandra 
observations set tight upper limits (99.7 per cent confidence level) in both the X-ray lumi- 
nosity of the posited AGN, L x (2 - lOkeV) < 3 x 10 40 ergs _1 , and the black hole mass, 
< 10 4 A -1 M Q , where A is the Eddington ratio. The X-ray luminosity upper limit is insuffi- 
cient to explain the high ionisation line [OIV] 25.89/im, which is observed in the mid-infrared 
spectrum of the MRK 996 and is proposed as evidence for AGN activity. This indicates that 
shocks associated with supernovae explosions and winds of young stars must be responsible 
for this line. It is also found that the properties of the diffuse X-ray emission of MRK 996 
are consistent with this scenario, thereby providing direct evidence for shocks that heat the 
galaxy's interstellar medium and contribute to its ionisation. 

Key words: galaxies: individual: Markarian 996 - galaxies: active - galaxies: dwarf - X-rays: 
galaxies 



1 INTRODUCTION 

Understanding how supermassive black holes (SBH) grow and how 
they relate to the formation of their host galaxies are challenges 
for modern astrophysics. Black holes with masses « 10 9 M have 
now been observed out to z ~ 6, when the Universe was only 
about 1 Gyr old (e.g. |Fan et aT]|2006[ >. Simulations can produce 
such monsters at early epochs by postulating violent high accretion 
rate events triggered by mergers of gas rich galaxies (e.g. Volontcri 
|& Rees||2005) . These simulations however, have to assume seed 
black holes, which form very early in the Universe (z « 20) and 
which will eventually grow to masses as high as 10 9 Mq via ac- 
cretion. The mass distribution of the seed black holes is an im- 
portant parameter of the simulations. The remnants of the first stars 
that formed in the Universe (Population I II) are believed to produce 



black holes with masses up to 10 Mq l Abel et 



al.|2000 
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|et al.|2002> . In this case however, accretion well in excess of the 
Eddington rate is required to grow those seeds to « 10 9 Mq in 
less than 1 Gyr (e.g. |Volonteri & Rees|2005||Begefman et al.|2006) . 
An alternative scenario for the formation of seed black holes is the 
direct collapse of pre-galactic gas discs (Lodato & Nataraja n|2006| 



2007 1, This model produces black holes with masses up to 10 M, 



orders of magnitude more massive than Population III stars, thereby 
significantly relaxing the requirements for super-Eddington accre- 
tion. A powerful probe of the seed black hole population at early 



epochs are intermediate mass black holes ( < 10 6 M©) in the local 
Universe, which are expected to reside in low mass galaxies (i.e. 
M BH - o relation; [Ferrarese & Merritt|2000||Gebhardt et al.|2000| 
|Greene et al.|2010[ l. |Volonteri et a l. (20081 for example, argue that 
the black hole occupation number of dwarf galaxies can constrain 
models of the origin of the seed black holes in the early Universe, 
e.g. Population III stars vs direct collapse of pre-galactic gas discs. 



The dynamical signatures of moderate size black holes are 
hard to detect with current instrumentation. As a result searches 
for such objects in nearby galaxies have turned to the easier but 
less direct approach of finding AGN in small or bulgeless galaxies, 
i.e. black holes during their active stage. NGC 4395 and POX 52 
are the best studied examples of nearby br oad-line AGN with esti 
mated black holes mass of about 10 Mq i Filippenko & Ho 



2003 



|Barth et aT 2004 Pet erson et al.|2003) . Furthermore, optical spec 
troscopic studies have recently found evidence for type-2 AGN ac- 
tivity associated with intermediate mass black holes in the late type 
spirals NGC 1042 jSeth et al.|2008{|Shields etaTp Q08 1, UGC 6192 
( |Barth et al.|2008l > and NGC 3621 ( |Barth et al.|2009] >. Systematic 
efforts to compile samples of moderate size black holes have also 
been carried out by |Greene & Hoi <2004| |2007[ > and |Dong et al.] 
( |2007^ . Those studies used the SDSS spectroscopic database to 
identify galaxies with AGN signatures in their optical spectra and 
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estimated black hole masses of < 10 6 M Q . The hosts of those AGN 
are typically late-type galaxies with stellar masses < 10 10 Mq. 

A special class of low mass galaxies which have attracted 
much attention over the years are Blue Compact Dwarfs (BCDs). 
These systems are exceptional because of their low metallicities 
(1/40-1/3 solar) and high specific star-formation rates, which in- 
dicate one of the first major episodes of star-formation. BCDs are 
therefore unique laboratories for testing theories of galaxy forma- 
tion. They resemble the first galaxies that formed in the early Uni- 
verse out of gas with near-primordial metallicity, and which accord- 
ing to the standard hierarchical formation paradigm are the build- 
ing blocks of massive galaxies. An important recent development 
in the study of BCDs has been the suggestion that they harbor AGN 
associated with intermediate mass black holes. If true, this would 
offer the opportunity to study in detail the parallel formation of the 
first stars and the seed black holes in systems that resemble proto- 
galaxies. 

Evidence for AGN activity in BCDs inlcude, the incidence of 
dense ionised cores (> 10 4 cm -3 ), br oad Balmer emissio n lines 
with widths in excess of 2000 km s -1 l lzotov et al.|2007 



Izotov 



& Thuan 2008 ), high ionisation narrow optical emission lines sug- 
gesting the presence of a hard ionising continuum ( |Thuan & Izotov| 
[20051 iJaines et al.|2009) . The evidence for AGN in BCDs however, 
is far from conclusive. Broad and high excitation lines can be due 
to stellar activity, e.g. shocks in fast winds of young massive stars 
(e.g. |Izotov et al.|2007[[Thuan et al.|2008{| James et al.|2009). Un- 



fortunately, optical emission line diagnostic diagrams (e.g. Baldwin 
et al. 1981 ; Kewley et al. 2001 1, which have traditionally been used 
to separate stellar photoionization-dominated from AGN-excited 
sources, break down at the low metallicity regime lGrove s~et al.| 
|2006[ >, and hence cannot be applied to BCDs. The most reliable 
method for confirming or refuting the presence of an AGN in any 
galaxy, irrespective of metallicity, are X-ray observations. The de- 
tection of an X-ray point source at the nuclear regions of a galaxy 
is strong evidence for accretion on a massive compact object. 

In this paper we explore the X-ray properties of MRK996, 
a BCD that is among the prime candidates for hosting an AGN. 
Mid-IR spectroscopy of MRK996 with Spitzer IRS shows the 
[OIV] 25.89 /im line with an ionisation potential of 54.9 eV, indi- 
cating the presence of a hard continuum ( T huan et al.|20 08 1, which 
could be produced by an AGN. Optical integral field spectroscopy 
by VLT/VIMOS have revealed an extremely dense nuclear region 
(electron density ~ 10 7 cm -3 ), much denser than in typical BCDs 
(James et al. 2009}. Also, a very red point source (J — K — 1.8) 
was found at the nucleus of the galaxy, while photometry with 
Spitzer indicates hot dust in the central galaxy regions and a mid- 
IR Spectral Energy Distribution (SED) similar to that of Seyferts 
(Thua n et aL]|2008) . New Chandra observations of MRK996 are 
presented in this paper and are used to search for AGN signatures in 
this galaxy. We adopt a distance of 22.3 Mpc for MRK 996, which 
corresponds to a velocity of 1642 kins -1 (James et al. 2009) and 
H = 73. 5 km s -1 Mpc -1 . 



2 DATA PROCESSING 

MRK 996 was observed for a total of 50 ks with the ACIS-S (Ad- 
vanced CCD Imaging Spectrometer) camera aboard Chandra in 
August 3, 2009. The ACIS-S was operated in very faint mode and 
the target was placed at the nominal aimpoint of the S3 CCD. The 
data were reduced using the Chandra Interactive Analysis of Ob- 
servations (CIAO) package version 4.1.2 and the methodology de- 



scribed in |Laird et al.| J2009) . The first step was to produce the 
level 2 event file from the raw data. Hot pixels and cosmic ray af- 
terglows were flagged. The charge-transfer inefficiency (CTI) and 
time-dependent gain corrections were then applied and the ACIS 
pixel randomization was removed. As MRK 996 was observed in 
very faint mode, the ACIS particle background cleaning algorithm 
was also applied. The data were inspected for high particle back- 
ground periods and it was found that they were not affected. Im- 
ages were constructed in four energy bands 0.5-7.0 keV (full), 
0.5-2.0 keV (soft), 2.0-7.0 keV (hard) and 4.0-7.0 keV (ultra-hard). 
Sources were detected using the methodology described in |Laird| 
|et al.|p&09) . The final source list consists of sources with Pois- 
son probabilities < 4 x 10 -6 in at least one of the energy bands 
defined above. Source free background maps and sensitivity maps 
were also constructed as described in Georgakakis et al. ( 2008 ). 

The astrometry of the Chandra data was refined by comparing 
the positions of X-ray sources with those of optical galaxies in the 
HST Wide Field Planetary Camera 2 (WFPC2) observations pre- 
sented by |Thuan et al.|(T996| >. The astrometric accuracy of the HST 
image is 0.25 arcsec. Excluding sources within the optical radius 
of MRK 996 (0.3 arcmin) obtained from the Third Reference Cata- 
logue of Bright Galaxies (RC3; [de Vaucouleurs et ahp9 91 ), there 
are only 2 X-ray sources within the HST/WFPC2 field of view. 
Both of them are associated with optical galaxies. Comparison of 
the optical and X-ray source positions shows small systematic off- 
sets of <5RA = 0.25 arcsec and <5DEC = 0.05 arcsec, which were 
applied to the X-ray source catalogue and images. 



3 X-RAY POINT SOURCES 

At the Poisson probability threshold of 4 x 10 -6 , two X-ray 
point sources are detected within the optical radius of MRK 996 
(0.3 arcmin). Figure 1 (left) overplots the positions of the two 
sources on the HST image of MRK 996. None of them is associ- 
ated with the nucleus of the galaxy. Relative to the optical centre of 
MRK 996, one of the point sources lies ~ 7 arcsec to the south and 
the other ~ 6 arcsec to the north. The southern X-ray detection is 
likely associated with one of the globular clusters of MRK 996. It 
lies m 0.15 arcsec from an optical point source on the HST image 
of the galaxy identified as a globular cluster by Thuan et al. (1996). 

The X-ray fluxes and luminosities of the two point sources 
are estimated from their X-ray spectra. The SPECEXTRACT task of 
CIAO was employed to extract the counts at the positions of each 
source, estimate background spectra, generate the Redistribution 
Matrix Function (RMF) and the Ancillary Response File (ARF). 
Source spectra were grouped to have at least one count per bin. The 
spectral analysis was carried out within the XSPEC vl2 package. 
As the net counts are low for both sources (see Table [TJ, we fit 
the data with a simple power-law model absorbed by cold gas in 
our Galaxy. A Galactic hydrogen column density of Nh = 4 x 



was adopted using the HI map of Kalberla et al. 



(20051 



In XSPEC terminology, our model is wabs x po, where wabs is 
the cold absorption from our Galaxy and po is the power-law. The 
derived parameters are listed in Table[T] The luminosities of the two 
sources are typical of X-ray binaries in galaxies. 

We assess the probability that the two off-nuclear X-ray 
sources are background AGN. We use the 0.5-10keV log JV— log S 
relation of Georgaka kis et ah] ( |2008jl to estimate the expected 
number of sources brighter than 10 -lb erg s -1 cm -2 within a ra- 
dius of ~ 7 arcsec from the optical centre of MRK 996. The flux 
limit above corresponds to the approximate 0.5-10keV flux of the 
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Figure 1. Left: HST F569W filter image of MRK 996. The circles are 0.5 arcsec in radius and mark the positions of the X-ray point sources detected in the 
Chandra observations. The contours correspond to the diffuse X-ray emission of MRK 996. North is up, East is left. Right: Zoom in the nuclear regions of 
MRK 996 to compare the spatial distributions of the hot X-ray emitting gas versus the optically detected HII gas. The greyscale is the HST F569W filter 
image, the red contours correspond to the diffuse X-ray component (same as contours on the left image) and the green contours are the narrow component of 
the Ho emission line from the integral field spectroscopic observations described by James et al. (2009). 



faintest of the two sources. The radius cut is the maximum angu- 
lar distance of the two point sources from the galaxy centre. The 
expected number of background sources is 0.12. The Poisson prob- 
ability of two background AGN within ~ 7 arcsec of the MRK 996 
nuclear regions is 1.8 x 1CP 4 . We conclude that both X-ray point 
sources are likely to be associated with MRK 996. 



4 DIFFUSE EMISSION 

The nuclear regions of the MRK 996 are dominated by diffuse X- 
ray emission. We explore the morphology of the X-ray gas by first 
removing point sources and then adaptively smoothing the full band 
X-ray image using the CMOOTH task of CIAO. The X-ray contours 
of the smoothed image are overplotted on the HST observations in 
Figure[T] The peak of the diffuse X-ray emission lies in between the 
nuclear star-forming region and a point source in the SW which was 
proposed to be a young star cluster by |Thuan et al.|(T9 96). There 
is also tentative evidence, limited by small photon statistics, that 
the distribution of the diffuse component is more extended in the 
southern direction. This is interesting, if it is confirmed by deeper 
X-ray observations, as the spatial distribution of globular clusters 
of MRK 996 also appears to be assymetric in the same direction. 
Most of them are found in the south of the galaxy iThua rTet al.| 
|1996| >. Figure [T] (right) also shows that the spatial distribution of 
the diffuse X-ray component and the H II ionised gas are similar. 

The X-ray spectrum of the diffuse component has been ex- 
tracted using SPECEXTRACT task of CIAO which is the proposed 
tool in the case of extended sources. The extraction aperture has 
a radius of about 2 arcsec and encompasses most of the diffuse 
component photons. The background spectrum was determined 
from nearby source-free regions. The spectrum was fit with a hot 
thermal plasma model Raymond & Smith] ( |1977) which is ab- 
sorbed by cold gas in our Galaxy with an HI column density of 
Nh = 4 x 10 20 cm -2 . The X-ray spectrum and the best-fit model 
are shown in Figure[2] Table[2]lists the hot gas parameters estimated 
from the best-fit model to the data. The temperature of the thermal 
plasma, about 1 keV (~ 10 7 K), is higher than the typical tempera- 



MRK 996 extended X-ray emission 
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Figure 2. Spectrum of the diffuse X-ray emission of MRK 996. The best 
fit thermal plasma model is shown with the histogram. The residuals of the 
model are also plotted. 



tures of the diffuse hot gas in dwarf starbursts (e.g. |Ott et al.|2005] > 
and other BCDs ( |Thuan et al.|2004> . 



5 DISCUSSION 

5.1 Is there an AGN in MRK 996? 

One of the motivations for the Chandra observations of MRK 996 
was to explore the possibility that an AGN is responsible for the ex- 
citation of the mid-IR [OIV] 25.89 /im line reported by Thu an et al.| 
(20081. Below about 2keV the nuclear regions of the galaxy are 
dominated by diffuse hot gas, which renders the search for a point 
source in the soft and full bands difficult. The hard band (2-7 keV) 
is less contaminated by the diffuse X-ray emission. We therefore 
estimate a conservative upper limit (99.7 per cent confidence level) 
for the flux of a point source at the optical centre of MRK 996 as- 
suming that the observed counts at that position are from a nuclear 
X-ray point source. 
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Table 1. MRK996 point source spectral parameters. Columns are (1) CXO ID, (2) total counts within the extraction region, (3) best-fit photon index of the 
power-law model, (4) c-statistic value of the fit and degrees of freedom, (5) flux in the 0.3-10keV band, (6) X-ray luminosity in the 0.3-10keV band. 



ID 


Total counts 


r 


c-stat/dof 


/ x (0.3-10keV) 
(lO-^ergs- 1 cm" 2 ) 


L x (0.3 - lOkeV) 
(ergs -1 ) 


CXOJ012735.6-061930 
CXOJ012735.5-061943 


74 
11 


4.0 ±0.6 

z -°-1.6 


33.7/24 
8.9/9 


2 0+ 10 
z - u -0.6 

3+ 02 


1.2+°-° x 10 39 

i.st ; 2 x io 38 



Table 2. The gas parameters for the thermal plasma model used to fit the 
spectrum of the diffuse emission of MRK996. The parameters are esti- 
mated for a metallicity Z opt = 0.5 Zq determined by James et al. 2009 
from the optical emission line spectrum of MRK996 and a volume fill- 
ing factor f v = 1. The parameter values can be converted to lower filling 
factors and other metallicities using the scaling relation given in the Ta- 
ble, where £ = Z/Z op t. The listed parameters are: (1) best plasma tem- 
perature in keV, (2) the normalisation of the best-fit model, which equals 
10~ 14 (4 7T D 2 )- 1 Jn e n H dV, where D is the distance to MRK996 in 
cm, n e , rij] are the electron and proton densities in cm -3 , (3) the plasma 
temperature in K, (4) The 0.3-10keV X-ray luminosity in erg/s, (5) the 
mean electron density in cm -3 estimated from the normalisation of the 
best-fit model assuming the proton and electron densities are equal, (6) pres- 
sure of the hot gas estimated from the relation P/k = 2n e T, (7) the gas 
thermal energy, E t h = 3 n e V kT, where V is the volume of the X-ray 
emitting region, which is assumed to be a sphere of 2 arcsec=214pc radius 
at the source distance (i.e. similar to the extraction region of the X-ray spec- 
trum of the diffuse component), (8) total mass of the hot gas estimated from 
the relation M^ot = Tie m p V, where m p is the proton mass, (9) gas cool- 
ing time, t coo i = Eth/Lx and (10) X-ray emissivity estimated from the 
relation Ax = Lx /(norm x4ttD 2 x 10 14 ). 



Parameter 


Value 


kT (keV) 


1 1±+ - 4 


norm (10 — 6 ) 


1 S+ ' 9 


T(10 7 K) 


1 o , +0.5 


L x (0.3 - 10 keV) (10 38 erg/s) 


1 1++ - 6 


<n e >[x(/0- OI3 ](cm- 3 ) 


o io +0 - 03 


P/k [x(/„£)- a5 ] (10 6 Kcm- 3 ) 


2 7+ 1 ' 4 
z - ' -0.7 


Eth [x/0- 5 £-°- 5 ](10 53 erg) 


tr n+3.6 
°- z -1.9 


M hot [x/»' 5 r°- 5 ](10 5 M Q ) 




tcooi [x/„ a5 r°- 5 ](10 8 yr) 


1 -°-0.4 


A (IO -23 ergs -1 cm 3 ) 


1 0+ ' 7 



For the estimation of the upper limits we adopt a Bayesian 
methodology similar to that described in |Laird et al.| ( |2009] >. For a 
source with observed number of total counts N (source and back- 
ground) within the 70 per cent EEF (Encircled Energy Fraction) 
radius and a local background value B the probability of flux fx is 

T N -T 

P(/x,JV) = — jvj— - (D 
where N = S + B and S = fx x t exp x C x q. In the last equation 



t exp is the exposure time at a particular position after accounting 
for instrumental effects (i.e. exposure map), C is the conversion 
factor from flux to count rate which depends on the adopted model 
for the spectrum of the source and 77 is the encircled energy fraction, 
i.e. 0.7 in our case. The integral of P(fx,N) gives the probability 
that the flux of the source is lower than fx,u, i.e. the confidence 
interval CL for the flux upper limit 

[ fXU P{fx,N)df x =CL. (2) 

J fx,L 

The lower limit of the integration is fx,L = 10 -18 erg s -1 cm -2 . 
Equation [2] is then solved numerically to estimate the flux upper 
limit that corresponds to the confidence interval of 99.7 per cent. 
For the flux to count rate conversion factor we adopt two extreme 
values. The first corresponds to a power-law with T = 1.9 and rep- 
resents the case of an unabsorbed AGN. The second is for a reflec- 
tion dominated spectrum, which is approximated by the PEXRAV 
model of XSPEC ( Magdziarz & Zdziarski 1995 ) and represents the 
case of a Compton Thick AGN (N H > 10 24 cm" 2 ). The 99.7 per 
cent upper limit is f x (2 - lOkeV) = 6 x 10" 15 ergs -1 cm -2 
for the r = 1.9 model and 10 -14 erg s _1 cm -2 for the reflection 
dominated spectrum. The 2-10 keV luminosity 99.7 per cent upper 
limits are respectively, 3.5 x 10 38 and 6 x 10 38 erg s _1 . In the case 
of a reflection dominated spectrum, we assume that the intrinsic lu- 
minosity of the source is 50 times higher (3 x 10 40 erg s _1 ), as the 
reflected emission is believed to represent 2 per cent of the direct 
component (e.g. Gilli et al. 2007). The corresponding upper limits 
for the black hole mass are < 100 A -1 M in the case of no ob- 
scuration and < 10 4 A -1 Mq for Compton thick obscuring clouds, 
where A is the Eddington ratio. In this calculation a bolometric to 
hard X-ray luminosity ratio of 35 is adopted (Elvi s et al.|1994| >. 

The luminosity upper limits above should be compared with 
the AGN luminosity required to explain the [OIV] 25.89 /im line in 
the mid-IR spectrum of MRK996 reported by |Thuan et~aT1 ( |2008) . 
Those authors modeled this line with an AGN ionization spectrum 
of the form f v oc v~ x and estimated a total rate of ionizing photons 
\ogQ(H) = 51.125. Assuming the photons are distributed in the 
energy interval 1-1000 Ry, this rate translates to a 2-10 keV intrin- 
sic AGN luminosity of 4.5 x 10 40 erg s _1 . This is 2 dex higher than 
the 99.7 luminosity upper limit in the case of an unabsorbed AGN. 
This indicates that an unabsorbed or mildly absorbed active SBH 
cannot be responsible for the ionisation of the [OIV] 25.89 ^tm line. 
If the central engine in MRK996 is absorbed by Compton Thick 
material and the X-ray emission is dominated by the reflected com- 
ponent, the estimated 99.7 per cent upper limit for the intrinsic lu- 
minosity of the AGN is still lower compared to that required to 
explain the [OIV] 25.89 fim line. We conclude that the X-ray data 
provide evidence against AGN ionisation for the [OIV] 25.89 71111 
line at a significance level of about 99.7 per cent. 
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5.2 Origin of the diffuse component 

An alternative explanation for the [OIV] 25.89 /an line is ionisation 
via shocks (Thuan et al. 2008 1 associated with stellar outflows and 
supernovae explosions. Shock excitation has already been proposed 
to partly explain the broad component of the optical emission lines 
in MRK 996 \ James et' al. 20091. The shocks generated by starburst 
winds are also known to carve bubbles of hot (10 6 - 10 7 K) and 
overpressured gas in the inter-stellar medium, which can be ob- 
served at soft X-rays (e.g. |Ott et al.|2005| >. 

We explore the possibility of a superbubble for the observed 
diffuse X-ray emission by comparing the observed X-ray luminos- 
ity of the hot gas, as determined from the X-ray spectral analysis of 
section [4| with the pre dictions of the bubble models of Castor et al. 
fl975) , |Weaver et al.|(T977) >, |Mac Low & McCray|(l988 i~~ 



n(r) 2 A x (T,Z)dV, 



(3) 



where Ax (T, Z) is the X-ray volume emissivity of the gas and 
n(r) is the radial density profile of the bubble. Substituting in 
the integral above the n(r) relation of Chu et al. 1 1995} and the 
Ax (T, Z) estimated from the observations (see Table |2|l we find 



2.6 x 10 34 JL 



-33/35 17/35,19/35 

n n tf- 



(erg s 



(4) 



where L37 is the mechanical luminosity of the starburst in units 
of 10 37 ergs _1 , no is the number density of the ambient ISM in 
cm -3 and ta is the age of the starburst in Myrs. The dimensionless 
integral / has a value of about 2 (see |Chu et aL]|1995| >. We do 
not include an explicit dependence on metallicity in equation [4] as 
in |Chu et al.| ( |1993] >. Such a dependence is accounted for in the 
calculation of Ax (T, Z) through the normalisation of the thermal 
plasma model fit to the spectrum of the diffuse X-ray component. 

For the estimation of the mechanical luminosity injected by 
the the winds of young stars and supernovae explosions we use 
the STARBURST99 stellar evolution synthesis model ( |Leithererj 
|et al.p 999; Vazquez & Leitherer 2005). We adopt a continuous 
star-formation model for the starburst with a Salpctcr IMF (lower 
mass IMq, upper mass 100 A/q, slope 2.35) and star-formation 
rate of IM0 yr -1 and metallicity Z = 0.008, i.e. similar to the 
0.5 Z Q metallicity of MRK 996 < |James et al.]|2009l >. The model 
mechanical luminosity is then scaled to the star-formation rate of 
MRK 996, for which we adopt the value 1.3 Mq yr - . This was es- 
timated from the integral field spectroscopy observations presented 
by I James et al.|p009| > by integrating the narrow line component 
of the Ha emission line under the assumption that the broad com- 
ponent of this line is associated with shocks. [James et aL| (2009) 
estimated an age for the young stellar population in MRK 996 of 
3-5 Myr based on the equivalent width of the H/3 line. An aver- 
age starburst age of 4 Myr is adopted in equation [4] Under those 
assumptions the mechanical luminosity from the STARBURST99 
model is 3.9 x 10 40 erg s _1 cm -2 (Figure 1 12 of Leitherer et al. 
1999). For the density of the ambient ISM James et al. (2009) esti- 
mate no = 1 — 40 cm" 3 for the outer regions of MRK 996 based 
on the [S II] A6716 / A6731 intensity ratio. In ionized nebulae how- 
ever, the [S II] emission is typically associated with lower ioniza- 
tion gas which would be expected to be of higher density than the 
average. It is therefore likely that the lower end of the no range 
determined by James et al. (2009) represents the mean ISM den- 
sity into which the superbubble expands. Adopting no = 1 cm -3 
we predict an X-ray luminosity of 3 x 10 38 ergs _1 , larger than 
the observed one, 1.1 x 10 38 erg s _1 . However, given the large un- 
certainties in the model parameters (e.g. a slightly younger age of 



3 Myr yields Lx = 7 x 10 37 ergs _1 ) and the observables (e.g. 
Ax(T, Z), metallicity of the X-ray emitting gas, ISM density) we 
conclude that the observed Lx of the diffuse component is broadly 
consistent with the expanding bubble model prediction. 

We caution that the diffuse component could be the result of 
the combined emission of high mass X-ray binaries (HMXRBs) 
that lie below the sensitivity limit of the Chandra observations. We 
assess the contribution of such sources by integrating the global lu- 
minosity function of HMXRBs in starburst galaxies determined by 
|Grimm et al. ( 2003 1. A star-formation rate of 1 .3 Mq / yr is adopted 
for MRK 996 j James et al.|2 009 1 and the integration is carried out 
between 10 37 and 10 a8 erg s _i . The upper limit is the approximate 
0.5-10keV band point source sensitivity limit of the Chandra obser- 
vations. We estimate a total X-ray luminosity for HMXRBs below 
the Chandra data detection threshold of ~ 10 39 erg s _1 . In princi- 
ple the diffuse X-ray component could be attributed to HMXRBs. 
The X-ray spectrum of the diffuse component is not inconsistent 
with this interpretation. The best-fit spectral index of the power-law 
component is 2.2 ± 0.7. Although HMXRBs typically have flatter 
X-ray spectra, ~ 1.2 (e.g. Whi te et al.|1983fr , systems with spec- 
tral indices as flat as F « 2 have been observed (e.g. |Sasaki et al.| 
|2003[ >. Although it is hard to reject the possibility of X-ray binaries 
dominating the diffuse component there are arguments against this 
interpretation. The age of the starburst in MRK 996 is 3 — 5 Myr 
l |James et al.||2009| >, while models of the evolution of the X-ray 
luminosity of HMXRBs (e.g. |Sipior|2003||Linden et al.|2010) pre- 
dict that these sources peak at later times, ~ 10 — 20 Myr after 
the initial burst, especially in the case of low metallicity systems. It 
is also interesting that the peak of the diffuse X-ray emission falls 
between the nuclear starburst and the young super-star cluster can- 
didate in the SW of the nucleus identified by Thuan et al. (1996). 
One could reasonably assume that the peak emission is associated 
with an interaction region where the respective stellar wind super 
bubbles from the two sources collide and heat the ISM. 



6 CONCLUSIONS 

Chandra X-ray observations are used to search for an AGN in 
the BCD MRK 996, which could explain the high ionisation line 
[OIV] 25.89 (im observed in the mid-IR spectrum of the galaxy. 
The estimated upper limit for the X-ray luminosity of a nuclear 
point source in MRK 996 is inconsistent with AGN ionisation being 
responsible for the excitation of the [OIV] 25.89 «m line. Shock ex- 
citation must be responsible for this line. We find direct evidence 
for shocks by studying the diffuse X-ray emission of MRK 996. 
The properties of this component are broadly consistent with shock 
heating from stellar winds and supernovae explosions. A tight up- 
per limit of < 10 4 A _1 Mq is also set for the posited black hole 
mass ofMRK996. 



7 ACKNOWLEDGMENTS 

AG acknowledges financial support from the Marie-Curie Reinte- 
gration Grant PERG03-GA-2008-230644. YGT acknowledges the 
award of a Marie Curie intra-European Fellowship within the 7 th 
European Community Framework Programme (grant agreement 
PIEF-GA-2009-236486). Support for this work was provided by 
the National Aeronautics and Space Administration through Chan- 
dra Award Number 1 1700503 issued by the Chandra X-ray Obser- 
vatory Center, which is operated by the Smithsonian Astrophysical 



6 Georgakakis et al. 



Observatory for and on behalf of the National Aeronautics Space 
Administration under contract NAS8-03060. 



REFERENCES 

Abel T., Bryan G. L., Norman M. L., 2000, ApJ, 540, 39 
Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP, 93, 5 
Barth A. J., Greene J. E., Ho L. C, 2008, AJ, 136, 1 179 
Barth A. J., Ho L. C, Rutledge R. E., Sargent W. L. W., 2004, 
ApJ, 607, 90 

Barth A. J., Strigari L. E., Bentz M. C, Greene J. E., Ho L. C, 

2009, ApJ, 690, 1031 

Begelman M. C, Volonteri M., Rees M. J., 2006, MNRAS, 370, 
289 

Bromm V., Coppi P. S., Larson R. B., 2002, ApJ, 564, 23 

Castor J., McCray R., Weaver R, 1975, ApJ, 200, L107 

Chu Y, Chang H., Su Y., Mac Low M., 1995, ApJ, 450, 157 

de Vaucouleurs G., de Vaucouleurs A., Corwin Jr. H. G, Buta 
R. J., Paturel G., Fouque P., 1991, Third Reference Catalogue 
of Bright Galaxies. Volume 1-3, XI, 2069 pp. 7 figs, Springer- 
Verlag Berlin Heidelberg New York. 

Dong X., Wang T., Yuan W., Zhou H., Shan H., Wang H., Lu 
H, Zhang K., 2007, in Astronomical Society of the Pacific Con- 
ference Series, Vol. 373, The Central Engine of Active Galactic 
Nuclei, L. C. Ho & J.-W. Wang, ed., pp. 57-+ 

Elvis M., et al., 1994, ApJS, 95, 1 

Fan X., et al., 2006, AJ, 131, 1203 

Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Filippenko A. V, Ho L. C, 2003, ApJ, 588, L13 

Gebhardt K., Bender R., Bower G., Dressier A., Faber S. M., Fil- 
ippenko A. V, Green R., Grillmair C, Ho L. C, Kormendy J., 
Lauer T. R., Magorrian J., Pinkney J., Richstone D., Tremaine 
S., 2000, ApJ, 539, L13 

Georgakakis A., Nandra K., Laird E. S., Aird J., Trichas M., 2008, 
MNRAS, 388, 1205 

Gilli R., Comastri A., Hasinger G, 2007, A&A, 463, 79 

Greene J. E., Ho L. C, 2004, ApJ, 610, 722 

— , 2007, ApJ, 670, 92 

Greene J. E., Peng C. Y, Kim M., Kuo C. Y, Braatz J. A., Impel- 
lizzeri C. M. V, Condon J. J., Lo K. Y, Henkel C, Reid M. J., 

2010, ArXiv e-prints 

Grimm H„ Gilfanov M., Sunyaev R., 2003, MNRAS, 339, 793 
Groves B. A., Heckman T. M., Kauffmann G., 2006, MNRAS, 
371, 1559 

Izotov Y. I., Thuan T. X., 2008, ApJ, 687, 133 

Izotov Y. I., Thuan T. X., Guseva N. G, 2007, ApJ, 671, 1297 

James B. L., Tsamis Y. G., Barlow M. J., Westmoquette M. S., 

Walsh J. R., Cuisinier F, Exter K. M., 2009, MNRAS, 398, 2 
Kalberla P. M. W., Burton W. B., Hartmann D., Arnal E. M., Ba- 

jaja E., Morras R., Poppel W. G. L., 2005, A&A, 440, 775 
Kewley L. J., Dopita M. A., Sutherland R. S., Heisler C. A., 

Trevena J., 2001, ApJ, 556, 121 
Laird E. S., et al., 2009, ApJS, 180, 102 

Leitherer C, Schaerer D., Goldader J. D., Gonzalez Delgado 
R. M., Robert C, Kune D. F, de Mello D. F, Devost D., Heck- 
man T. M., 1999, ApJS, 123, 3 

Linden T, Kalogera V, Sepinsky J. F, Prestwich A., Zezas A., 
Gallagher J., 2010, arXivl005.1639 

Lodato G, Natarajan P., 2006, MNRAS, 371, 1813 

— , 2007, MNRAS, 377, L64 

Mac Low M., McCray R., 1988, ApJ, 324, 776 



Magdziarz P., Zdziarski A. A., 1995, MNRAS, 273, 837 
Ott J., Walter F, Brinks E., 2005, MNRAS, 358, 1453 
Peterson B. M., Bentz M. C, Desroches L., Filippenko A. V, Ho 
L. C, Kaspi S., Laor A., Maoz D., Moran E. C, Pogge R. W, 
Quillen A. C, 2005, ApJ, 632, 799 
Raymond J. C, Smith B. W., 1977, ApJS, 35, 419 
Sasaki M., Pietsch W, Haberl F, 2003, A&A, 403, 901 
Seth A., Agueros M., Lee D., Basu-Zych A., 2008, ApJ, 678, 1 16 
Shields J. C, Walcher C. J., Boker T, Ho L. C, Rix H., van der 

Marel R. P., 2008, ApJ, 682, 104 
Sipior M. S., 2003, PhD thesis, "Population synthesis and its con- 
nection to astronomical observables", The Pennsylvania State 
University 

Thuan T. X., Bauer F. E., Papaderos P., Izotov Y. I., 2004, ApJ, 
606, 213 

Thuan T. X., Hunt L. K., Izotov Y. I., 2008, ApJ, 689, 897 
Thuan T. X., Izotov Y. I., 2005, ApJS, 161, 240 
Thuan T. X., Izotov Y. I., Lipovetsky V. A., 1996, ApJ, 463, 120 
Vazquez G. A., Leitherer C, 2005, ApJ, 621, 695 
Volonteri M., Lodato G., Natarajan P., 2008, MNRAS, 383, 1079 
Volonteri M., Rees M. J., 2005, ApJ, 633, 624 
Weaver R, McCray R., Castor J., Shapiro P., Moore R., 1977, 
ApJ, 218, 377 

White N. E., Swank J. H., Holt S. S., 1983, ApJ, 270, 711 



